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Starting from the infinite-dimensional Ikeda map, we derive an extended temporal Lugiato-Lefever 
equation that may account for the effects of the conjugate electromagnetic fields (also called ‘negative 
frequency fields’). In the presence of nonlinearity in a ring cavity, these fields lead to new forms of 
modulational instability and resonant radiations. Numerical simulations based on the new extended 
Lugiato-Lefever model show that the negative-frequency resonant radiations emitted by ultrashort 
cavity solitons can impact Kerr frequency comb formation in externally pumped temporal optical 
cavities of small size. Our theory is very general, is not based on the slowly-varying envelope 
approximation, and the predictions are relevant to all kinds of resonators, such as fiber loops, 
microrings and microtoroids. 


PACS numbers: 42.65.-k, 42.65.Tg, 42.60.Da, 42.65.Sf 

Introduction — Kerr frequency combs (KFCs), i.e. 
light sources with a large number of highly resolved and 
nearly equidistant spectral lines, are attracting a consid¬ 
erable interest in recent years, due to their important 
applications in metrology, optical clocks, precision spec¬ 
troscopy, precision time and distance measurements, at- 
tosecond pulse generation and complex nonlinear dynam¬ 
ics, to name just a few m- Physical systems that are 
able to generate KFCs are microring resonators, micro¬ 
toroids, crystalline resonators, microspheres, photonic- 
crystal cavities and optical fiber loops [1]. Microring 
resonators are particularly important in this respect, 
since they are small-size, low-loss, CMOS-compatible and 
power efficient devices that can be made of different non¬ 
linear materials and are therefore ideal for on-chip KFC 
generation. 

The main ingredients for an efficient KFC formation 
have been identified to be four-wave mixing (FWM) and 
temporal cavity soliton (CS) generation [2]. There is cur¬ 
rently an intense research activity aiming to maximise the 
spectral extent of the comb and its coherence, and to un¬ 
derstand the experimentally obtained spectra from first 
principles. Due to the extremely complex dynamical be¬ 
haviour and stability properties of the propagating CSs 
and patterns in the resonators, an intense theoretical ac¬ 
tivity on the mathematical properties of the traditionally 
used averaged propagation equation, called the temporal 
Lugiato-Lefever equation (LLE), has developed over the 
past years, with a frequent display of new and surprising 
results [Sd- 

One of the major dynamical effects in the propagation 
of ultrashort pulses is the radiation emitted by solitons 
due to higher-order dispersion effects, also called reso¬ 
nant radiation (RR) jSUTT]. This radiation, which ap¬ 
pears when pumping near the zero-dispersion point of 
the structure, is very visible in experiments performed 
with optical fibers and it also plays a central role in the 
dynamics of CSs [SllIHil]. Indeed, it has been shown 
experimentally that CSs emit RR in microring resonators 


and fiber loops mns]. 

In addition to the traditional RR emission, other emis¬ 
sions are possible when considering the influence of conju¬ 
gate fields and the nonlinear interaction between positive 
and ‘negative frequencies’ in the nonlinear polarization 
HU. Recently, experiments performed in optical fibers 
and bulk crystals revealed the presence of ‘negative fre¬ 
quency resonant radiation’ (NRR), which has also been 
completely explained theoretically in Ref. [18] . The the¬ 
oretical approach used was one based on the analytic sig¬ 
nal, which is able to take into account the dynamics of 
both positive and negative frequencies, thus avoiding the 
use of the slowly-varying envelope approximation which 
breaks down for very short pulses [T8H2Q] . An intriguing 
question is whether nonlinear optical cavities may some¬ 
how enhance, and therefore be affected by NRR-like ef¬ 
fects. 

In this paper we extend the temporal LLE in order 
to take into account the effects of negative frequencies 
and conjugate fields on the propagation of CSs in opti¬ 
cal resonators {extended LLE, eLLE for brevity). This 
extension is necessary in order to be able to study ana¬ 
lytically and numerically the existence of new resonant 
radiations emitted by the CSs during the circulation in 
the resonator. We show that due to the forced and dissi¬ 
pative nature of the eLLE, the new resonant radiations, 
that are typically remarkably feeble in conventional op¬ 
tical fibers and bulk materials, can become quite strong 
and can be more efficiently generated. This surprising 
result has the potential to impact considerably the for¬ 
mation of KECs due to CSs and also the intrinsic stabil¬ 
ity of the homogeneous steady state CW solutions of the 
cavity. 

Extended temporal Lugiato-Lefever equation — The 
starting point of any discussion on optical resonators and 
cavities is the infinite-dimensional Ikeda map Eung 
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FIG. 1: (Color online) Sketch of the generic microcavity and 
its parameters, r is the ring radius, Pin is the cw input power, 
Pout is the output field power, is the intracavity field 

at position ^ in the ring, R and T are the reflection and trans¬ 
mission coefficients, respectively. 


Here, Ai^ is the envelope of the (impulsed or cw) pump 
field, T is the transmission coefficient at the coupling 
point z = 0 , P is the reflection coefficient at the same 
point (P^ + = 1 ), 00 = ^0 “ 27 rm is the phase accu¬ 

mulated over a round-trip, ( 5 o is the cavity detuning, L is 
the length of the cavity, An{z^ t) is the envelope of the in¬ 
tracavity field circulating at the n-th step, 2; is the spatial 
coordinate along the cavity, t is the ‘fast’ time variable 
in the reference frame moving at the group velocity 1/01, 
D{idt) = i^ dispersion operator, 

Pj = [dlj^{uj)]uj=ujo is the j-th dispersion coefficient cal¬ 
culated at the pump frequency cjo, S{idt) = (1 -{-iujQ^dt) 
is the shock operator describing the first order correc¬ 
tion due to the frequency dependent nonlinearity. For 
simplicity, we do not consider here the influence of the 
Raman effect, as this can also be easily included in our 
model. 

The Ikeda map 0i is usually based on the slowly- 
varying envelope approximation (SVEA), under which 
the spectral extent of the pulse must be much smaller 
than its central frequency. In Eq. ([^, Pni is the non¬ 
linear polarization of the intracavity field, which in pres¬ 
ence of the Kerr effect only can be written as Pn\[An] = 
j\An\‘^An, where 7 is the nonlinear coefficient of the ma¬ 
terial. However, as was explained in Ref. m, this simple 
form of the polarization does not describe sufficiently well 
the dynamics of ultrashort pulses. In order to include in a 
physically consistent way the higher-order effects such as 
third-harmonic generation (THG) and the contribution 
of conjugate terms, one must use the so-called analytic 
signal (AS) for all the fields involved in Eqs. 01)- We 
refer the reader to Refs. (TSHSO] foi* a full discussion, but 
we review here the salient features of ASs for complete¬ 
ness. 

The full real electric field propagating in an optical 
system is denoted by E(z,t), where 2; is the propagation 
direction and t is the time variable. Its Eourier transform 
is denoted by E^^{z) = E[E{zA)] = E{zA)E^^dt. 
The analytic signal of the electric field, i.e. the posi¬ 


tive frequency part of the field, is defined as S{zA) = 
E^{z)e~'^^^duj, and is a complex function. The 
Eourier transform of the electric field can be written as 
the sum + (f_^)*]/2 since only the positive (or 

negative) frequency part of the spectrum carries infor¬ 
mation, while for the same reason the electric field itself 
is real and is given by E{zA) = [£{zA) + 

The ‘envelope’ we introduce here is defined by A{z, t) = 
f (^, (where 0o is the propagation constant 

at the pump frequency), i.e. the frequency components 
of the AS are ‘shifted’ by an amount —cjo- By doing this, 
we deal with frequency detunings Acj from cjo, and not 
with absolute frequencies. However, there is a key differ¬ 
ence between the conventional definition of envelope and 
the envelope of the AS: the former is adequate only if the 
spectral extension of the pulse evolution is much smaller 
than the pulse central frequency, \Auj\ = |u; — cjol ^ cjo, 
i.e. only under SVEA conditions, while the envelope 
of the AS A{zA) considered here does not suffer from 
this limitation. By clearly dividing the envelope asso¬ 
ciated to the positive frequency components from that 
associated to the negative ones, we will be able to write 
the extended temporal LEE that correctly describes the 
dynamics of intracavity pulses of arbitrary duration and 
spectral extension^ taking into account the peculiar and 
non-trivial interaction between positive and negative fre¬ 
quencies that arises in the complete nonlinear polariza¬ 
tion. In this paper, all fields, linear and nonlinear, are 
intended to be envelopes of ASs. 

Within the above framework, the AS of the nonlinear 
polarization is written as m- 


Pn\ [^n] — T 




( 3 ) 

where (f^zA) = + Akz^ and Ak = {/3iUJo — 0 o) is a 

factor, crucial for the efficient phase-matching of the res¬ 
onant negative-frequency terms, that measures the dif¬ 
ference between the phase and the group velocity in the 
medium. The subscript + in Eq. © signifies the filtering 
of the negative frequency components out of the polar¬ 
ization. This ensures that during its evolution An{zA) 
only contains positive frequencies, and it is thus consis¬ 
tent with its own definition. The first term in Eq. © is 
the usual Kerr term, the second is the so-called ’negative 
Kerr’ (NK) term, while the third gives THG. Note that 
NK and THG must appear together in order to preserve 
the Hamiltonian nature of the EWM interaction, some¬ 
thing that was overlooked prior to our theoretical work 
on the subject [ 18 ]. 

The accurate numerical integration of Eqs. ( [Tp] ) is 
very time consuming, and the equations are not suitable 
for a deep analytical understanding. This is the reason 
why mean-field models based on averaging over many 
roundtrips (see Ref. [ 23 |) are a common tool for dramat¬ 
ically reducing the complexity of the map into a single 
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equation. To this aim, we now plug Eq. ^ into Eq. 
® , and we follow the averaging procedure described in 
[2^ . using the ‘cavity soliton units’ ^ = z/I/D 2 , r = t/tQ^ 
Q = iuoto, K = A/cI/d 25 ^0 = (7^02)”^ 

(so that Lnl = [ 7 ^ 0 ]“^ = ^ 02 ), V^in = An/v^, 
r) = Clr + A natural value for the time scale 
to is the typical single-CS duration to = [\f32\L/{26o)]^^‘^. 
In the averaging procedure, the cavity length L must be 
much smaller than the dispersive and nonlinear lengths, 
I/D 2 and I/nl respectively. This ensures that the in¬ 
tracavity pulse does not change much during a single 


roundtrip. This slow variation must also be satisfied by 
the NK and THG terms in Eq. (§. This means that 
the phase (t){z^t) must rotate very rapidly in a single 
roundtrip, so that the average effect is mediated almost 
to zero, resulting in the condition k ^ meaning that 
the group velocity and the phase velocity must be suffi¬ 
ciently different, which is true in dispersive media such 
as dielectric waveguides. 

In this way, we arrive at the following extended 
Lugiato-Lefever equation (eLLE): 


+ D{idr)i^ + '^(r + id)^lj + (1 + ^dr) 




- ifiXpin = 0 . 


(4) 


-I + 


T = [{a,L + T 2 )/ 2 ]Ld 2 /A, ^ = ^ 0 ^ 02 /^, /i = 

TLd 2 /A, D{idr) = Y.n> 2 ^n{idrY, where bn = 
/ 3 ^/(n!tQ“^|/ 32 |) are the dimensionless dispersion coeffi¬ 
cients, and r) = Or + 

Equation @ is the central result of this paper. It can 
be applied to any kind of optical resonator, and can be 
easily extended to include the Raman effect or any other 
perturbation of the NLSE. 

Phase matehing eonditions for the new radiations — 

We now derive the phase-matching conditions for the 
most important resonant radiations emitted by the ultra- 
short CSs propagating in the resonator. In Eq. Q, we 
substitute the Ansatz '0(r, t) = fjoP r), and linearize 
with respect to the small radiation field g. Note that due 
to the loss term proportional to T in Eq. Q, far from 
the central peak of the CS the resonant radiations decay 
asymptotically towards the complex cw background fjQ. 
Thus the intensity of the CS peak power does not play 
any important role, and the problem becomes very sim¬ 
ilar to a modulational instability (MI) phase-matching. 
However, the usual MI analysis does not work here, since 
the presence of the exponential terms in Eq. Q cannot 
give a stationary state, and the standard techniques used 
in Refs. [T2l|2l]|22] would not be appropriate. We there¬ 
fore use the perturbation method employed in Ref. [18] 
in the fiber-optics/bulk context, obtaining 

D{A)-vA = Do, (5) 

D{A) -vA = Do± 2n, ( 6 ) 

where A is the dimensionless frequency detuning from 
the pump, D{A) = X]n >2 is the dispersion of the 

linear waves, Dq = (5—iT— 2 |' 0 op is the complex nonlinear 
wavenumber, and v is the velocity parameter of the re¬ 
sulting moving CS, which starts to drift as a consequence 
of the higher-order dispersion m- Equation ^ is con¬ 
nected to the emission of the conventional RR, routinely 


observed in fibers and cavities nni HU E]. Equations 
phase-match the negative frequeney resonant radia¬ 
tion [NRR, ‘+’ sign in Eq. ®] , which is due to the NK 
term in Eq. © , and the third-harmonie resonant radi¬ 
ation [THRR, sign in Eq. ©] , which is due to the 
THG term in Eq. ©• 

Note that in all physically relevant situations k ^ 
The roots of Eqs. have real and imag¬ 

inary parts, indicating the frequency detuning of the 
emissions and their decay rate towards the soliton back¬ 
ground, respectively. Equation © is valid only in the 
presence of a propagating CSs, while Eqs. ^ holds 
also for a pure cw intracavity field. This latter phase¬ 
matching equation implies that the ‘old’ stable branch 
of the stationary solutions of Eq. 0, i.e. the smallest 
root of the equation i(r P i5)f; P |' 0 r '0 — = 0 (see 

also Refs. mm), is unstable in presence of NK and 
THG, and sidebands will appear. This effect, which is 
not present in the traditional temporal LEE, is also a 
new feature of our eLLE model, as we shall see shortly. 

Numerieal simulations — In order to illustrate the 
spectral dynamics and the emission of the resonant radi¬ 
ations from ultrashort CSs as seen in the previous section, 
we use a highly-nonlinear, small, low-loss resonator (for 
example a microtoroid or a microring, of the kind de¬ 
scribed in Refs. mm see also Eig. Q. We take a radius 
r = 30 /im, cavity length L = 1.88 x 10“^ m, pump wave¬ 
length Ao = 1.55 /im, P 2 = —90 psec^/km, f3s = — 1.11 
psec^/km, nonlinear coefficient 7 = 1 W“^m“^, group 
index at the pump Ug = 1.5, cw pump power Pin = 265 
mW, transmission coefficient T = 0.07, detuning from 
the cavity resonance (5o = 0.0115<C7r, free spectral range 
ESR=1060 GHz, roundtrip time tR = 1/ESR= 150 fsec, 
photon lifetime tph — 0.19 nsec, finesse P 628 and 
loaded Q-factor Q ^ 10^. The typical temporal width of 
the CSs formed in the cavity is to = [\f32\L/{26o)]^^‘^ — 27 
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fsec, and we take this value for the scaling of Eq. The 

typical peak power of the CS is also given by the soliton 
power scale Pq = 2(5o(I/D2/P) = [ 7 ^ 1 ) 2 ]”^ — 122 W. The 
second order dispersion length is I/d 2 = 8.2 mm, which 
gives a ratio I/D 2 /P — 43.5 ^ 1, making the averaging 
procedure meaningful. We take a value Ak ^ 0.67 x 10^ 
m“^, an order of magnitude that is common for solid 
media. It is very favourable to use small-size microcavi¬ 
ties in order to minimize the value of 0 accumulated over 
a roundtrip, proportional to A/cI/d 2 * ( 4 ^/ 4 ^D 2 ) = AkL. 
This gives the opportunity to observe the new nonlinear 
effects in realistic systems. 

With these parameters, the dimensionless coefficients 
in Eq. Q become T = 0.217, 6 = 0.5, /i = 3, 'ipin = 0.047, 
b 2 = ±0.5, 63 = -0.0758, O = 33 and - 54.5. Due 
to the fast oscillations in the exponentials in Eq. 
which can generate high-frequency radiations, we make 
sure to use a large number of sampling points in the split- 
step Eourier solver {N = 2^^), and a small spatial step 
(A^ = 10-^). 

Eigure [^a) shows the output spectrum of the propa¬ 
gation = 27, equivalent to 1200 cavity roundtrips) of 
the lower-branch homogeneous steady state (HSS) cw so¬ 
lution (see Refs. [24l[25]) in both anomalous and normal 
dispersions respectively (62 = ±0.5), by using the eLLE, 
Eq. 0 - when 63 = 0 . As discussed in the previous sec¬ 
tion, for these parameters the HSS solution is stable in 
absence of NK and THG terms. Conversely, accounting 
the complex conjugate fields leads to far-detuned side¬ 
bands (indicated by PI and P 2 in the figure) accord¬ 
ing to the phase-matching conditions of Eq. 0 . Note 
that this new form of instability appears irrespective of 
the sign of 62 , due to the ± in Eq. Moreover, the 

sidebands would be perfectly symmetric with respect to 
the pump frequency for vanishing higher-order dispersion 
terms, but become asymmetric when the contribution of 
bn >3 is important, see Eig. [^b). In this case, the two 
peaks are always slightly imbalanced in amplitude, due 
to the presence of THG and pump and loss in the sys¬ 
tem. This proves, for the first time to our knowledge, 
the existence of novel kinds of modulational instability 
that are due to the contribution of negative frequencies 
in the nonlinear polarization [Eq. (§]■ This MI can be of 
the order of 1 % of the input pump or larger and should 
therefore be readily visible in experiments. 

Eigure [^a) shows the formation of a stable moving 
CS in the cavity, in presence of third-order dispersion, 
63 7 ^ 0 and in anomalous dispersion 62 = —0.5, when 
using the eLLE Eq. 0 . The input seed is taken to be 

' 0 (^ = 0 ,r) = v^sech i.e. the autosolution of 

the pulsed cavity [25], which is a good approximation (al¬ 
beit with vanishing background) of the final CS. Eigure 
[^b) shows the real part of the phase-matching curves 
of Eqs. ([5][^, and the prediction of the frequency posi¬ 
tions of all the resonant radiations. Green dashed line 



FIG. 2: (Color online) (a) Output spectrum of the propaga¬ 
tion of the HSS cw solution for normal (^2 = +0.5, blue solid 
line) and anomalous (^2 = —0.5, green dashed line) disper¬ 
sions, in the case bs = 0. Propagation distance is ^ = 27. 
(b) Same as (a) but with bs = —0.0758. All other relevant 
parameters are given in the text. 

in Eig. [^c) shows the intracavity field spectrum, in the 
case when only the Kerr effect is present in the nonlinear 
polarization of Eq. 0 (i .e. the case of the conventional 
temporal LLE). In this case, only the usual RR peak, 
which satisfies Eq. 0 , is emitted near the pump (with 
a very small contribution of its symmetric counterpart, 
not discussed here, see [l 2 ]). Blue solid line in Eig. [^c) 
shows the intracavity field spectrum when taking into ac¬ 
count all the terms in the full polarization, which contains 
also the NK and the THG terms, see Eq. Q. One can 
notice the appearance of relatively strong peaks, which 
are the NRR and the THRR peaks described in the pre¬ 
vious section, and that satisfy the phase-matching condi¬ 
tions Eqs. 0 . Complete agreement is observed between 
the predicted radiation positions [Eig. [^b)] and the sim¬ 
ulated ones. Also note the appearance of the THG peak, 
located at 1 + A/D = 3. Eigure |^d) shows the position 
of the NRR and THRR peaks for different values of k. 
The amplitudes of the radiations are uneven since the CS 
background oscillates during propagation. 

Discussion and conclusions — Our results in this 
work show that frequency comb formation in Kerr me¬ 
dia is affected by the resonant radiations and MI result¬ 
ing from the nonlinear interactions between positive and 
negative frequencies. These new types of radiations play 
a conceptually important but somewhat minor role in 
optical fibers and bulk materials. However, in optical 
resonators such as fiber loops and microrings, the field 
is forced to circulate many times in the cavity, and CSs 
possess a cw background that may dramatically enhance 
the emission of the new radiations. As a consequence, 
the output KEG spectra are affected by the radiation 
peaks emitted by CSs, and the ‘negative frequency’ ra¬ 
diations play an important role. We have provided the 
derivation of a universal model (the eLLE), not based on 
SVEA, that is able to take into account the full dynamics 
of the ultrashort intracavity pulses, and is amenable to 
analytical investigations. 
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FIG. 3: (Color online) (a) Formation of a stable moving CS 
in the cavity, after a propagation of ^ = 27. (b) Rea l part of 
phase-matching curves for the three equations ( |5|6| ): (i) for 
Eq. (ii) for Eq. with + sign, and (iii) for Eq. with 
— sign. Circles are the predicted positions of all the generated 
radiations (RR, NRR, THRR). (c) Intracavity held spectrum 
at ^ = 27, when considering only the Kerr term (green dashed 
line), and when including the NK and the THG terms in Eq. 

(blue solid line). The resonant radiations generated by the 
NK and the THG terms are indicated with NRR and THRR, 
respectively, (d) Zoom of the NRR and THRR peaks when 
K = 54.5,65,76,87 (indicated by 1,2,3,4 respectively). The 
parameters of the simulation can be found in the text, and 
b2 = -0.5, bs = -0.0758, = 0.084. 


Different alternative (but equivalent) approaches ex¬ 
ist in the description of KFC formation [26l [27] . In one 
scheme, the dynamics of the annihilation operators of the 
cavity modes is described via Heisenberg equations. In 
the same spirit of our formulation, the inclusion of the 
negative-frequency effects could be accounted for by us¬ 
ing the complete (and correct) interaction Hamiltonian 
V = —{gh)[e^e^ee/2 + (e'^'eee + e’*'e'*'e’*'e)/3], in the lan¬ 
guage of Ref. [26|. The first term is the usual FWM 
scattering, while the second and third terms, which must 
come together since they are individually non-Hermitian, 
account for THG and NK effects, respectively. It would 
be interesting to see how the different approaches in the 


literature may converge towards a unified description of 
the phenomena described here. 


F.B. would like to thank S. Trillo, N. Joly, M. Erkin- 
talo, M. Clerici and L. Caspani for useful discussions. 
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